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Abstract 

In this paper we show that with standard methods it is possible to obtain highly precise results 
for QNMs. In particular, secondary modes are obtained by numerical integration. We compare 
several results making a detailed analysis. 
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I. INTRODUCTION 



Wave propagation around black holes is an active field of research (see The per- 
spective of gravitational waves detection in a near future and the great development of 
numerical general relativity have increased even further the activity on this field 2[]. Grav- 
itational waves should be especially strong when emitted by black holes. The study of the 
propagation of perturbations around them is, hence, essential to provide templates for the 
gravitational waves identification. Thus, activity in this field is developping quickly [3|. 

The Schwarzschild black hole is rather well known [l, Q] but it is important to get re- 
assured about the robustness of the methods and results. Therefore, here, we embark in a 
detailed study of the secondary modes by means of the subtraction of the first modes in time 
domain. The method, though very simple has not been used in the literature due to numeri- 
cal errors implicit in the method. However, we were able to use simple standard methods to 
show that the results are valid implying an acuracy of seven figures for the dominant mode 
and sometimes three for the secondary one. 

We have used the geometric system of units, for which h = c = G = 1. This means that 
the masses have dimension of length and are measured in metres. The conversion factor 
from metres to kilograms is c^/G. 



II. FUNDAMENTAL MODES AND FIRST OVERTONES 

We begin with a series of tables on the frequencies of the QNMs for the Schwarzschild 
black holes. We have employed, throughout our discussion, M = 1.0m since the black-hole 
frequencies scales as Muj = const. For the grid spacing, we have used h = lO^^m. In what 
follows, we have used the notation uj^^^'^ for the dominant mode, uj^^*-^ for the first overtone 
and, whenever applicable, uo'^^^- for the second overtone. 

A few remarks are due: the scalar £ = field oscillated almost nothing, so no frequencies 
were computed for it at all. The second overtone wasn't clearly visible for £ < 6, and even 
so only for £ = 8 and higher we could make minimally acceptable fittings to it. 

The figures ([1]) and ([2]) show that the second overtone appears clearly only for very high 
£ (£ > 6), and it is the reason why we have blank spaces for this overtone in tables ^ and 
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^DOM 


Si^DOM 


i^SEC 


^TER 


1 










2 


0.483644 - 0.0967590i 


3E-7 


0.470 - 0.300? 




3 


0.675367 - 0.0964997i 


2E-7 


0.667 - 0.288? 




4 


0.867417 - 0.0963923Z 


2E-7 


0.859 - 0.287? 




5 


1.059614 - 0.0966337i 


6E-7 


1.050 - 0.283i 




6 


1.251891 - 0.0963060Z 


9E-8 


1.224 - 0.283? 




7 


1.444214 - 0.0962866? 


7E-7 


1.433 - 0.282? 




8 


1.636565 - 0.0962724Z 


IE -6 


1.629 - 0.288? 


1.584 - 0.422Z 


9 


1.828959 - 0.0962639i 


IE -6 


1.822 - 0.288i 


1.791 - 0.404i 


10 


2.021329 - 0.0962568i 


IE -6 


2.014 - 0.288? 


1.990 - 0.405i 


11 


2.213730 - 0.0962522i 


IE -6 


2.207 - 0.288? 


2.174 - 0.440Z 


12 


2.406139 - 0.0962487Z 


IE -6 


2.400 - 0.288? 


2.371 - 0.440Z 



TABLE I: Frequencies for the Schwarzschild BH of M = 1.0m. Scalar Field, different i values. 





^DOM 




i^SEC 


i^TER 


1 


0.248229 - 0.0924905? 








2 


0.457595 - 0.0950044? 


hE-8 


0.440 - 0.290? 




3 


0.656899 - 0.0956165Z 


2E-8 


0.648 - 0.286? 




4 


0.853096 - 0.0958605i 


2E-7 


0.844 - 0.285i 




5 


1.047915 - 0.095982H 


9E-8 


1.039 - 0.282? 




6 


1.242000 - 0.0960523? 


9E-8 


1.232 - 0.282? 




7 


1.435647 - 0.0960959Z 


7E-7 


1.424 - 0.281? 




8 


1.629012 - 0.0961250i 


QE-8 


1.621 - 0.288i 


1.577 - 0.425i 


9 


1.822180 - 0.0961452i 


lE-% 


1.815 - 0.288? 


1.788 - 0.407i 


10 


2.015214 - 0.0961596Z 


lE-Q 


2.008 - 0.288? 


1.984 - 0.410Z 


11 


2.208148 - 0.0961714Z 


2E-6 


2.199 - 0.285? 


2.164 - QAA2i 


12 


2.401004 - 0.0961800i 


IE -6 


2.395 - 0.288? 


2.365 - 0.427Z 



TABLE II: Frequencies for the Schwarzschild BH of M = 1.0m. Electromagnetic Field, different £ 
values. 
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i 


^DOM 




^SEC 


^TER 


2 


0.37367 - 0.08896i 


6^-6 


0.352 -0.271i 




3 


0.599444 - 0.0927031i 


5^-9 


0.586 - 0.278i 




4 


0.809180 - 0.0941643i 


AE-?> 


0.797 - 0.279i 




5 


1.012297 - 0.0948713i 


3E-7 


1.002 - 0.279i 




6 


1.212013 - 0.0952667i 


2E-7 


1.203 -0.286i 




7 


1.409741 - 0.0955106i 


2E-7 


1.401 - 0.286i 


1.361 - 0.432i 


8 


1.606202 - 0.0956724i 


3E-7 


1.594 - 0.280i 


1.474 - 0.515i 


9 


1.801801 - 0.0957836i 


IE -6 


1.795 -0.287i 


1.770 - 0.41H 


10 


1.996796 - 0.0958649i 


2E-7 


1.990 -0.286i 


1.967 - 0.408i 


11 


2.191345 - 0.0959263i 


2E-7 


2.188 -0.287i 


2.154 - 0.415i 


12 


2.385555 - 0.0959724i 


2E-7 


2.379 - 0.287i 


2.353 - 0.446i 



TABLE III: Frequencies for the Schwarzschild BH of mass M = 1.0m. Axial Field, different i 
values. 

We have compared our data to those generated via 6th-order WKB computations, from 
4] and 3rd-order WKB from j^. When it comes to notation, we have employed £ for the 
multipole index (as before), n for the overtone index (n = for the fundamental, n = 1 
for the first overtone, and so on), ujnum refers to our numerical data, ujwkbi refers to the 
WKB data from paper [4] and ujwkb2, to paper j^. 

These comparison tables leave no room for doubts when it comes to the fundamental 
mode, n = 0. For all fields and ^-values, the data from our numerical simulations and 
those from the 6th-order WKB method showed a high-degree agreement, with differences 
between 1 part in 10^ and 1 part in 10^, for both Re{uj) and Im{uj), especially for higher 
^-values. When it comes to the first overtone, for all fields under study, the agreement 
between both sets of data was not so impressive, hovering around 1 — 2% for Re{uj) and 
2 — 3% for Im{ijj) for low ^ (typically up to £ = 6) and improving for higher £, to a few parts 
in 1000 for Re{uj) and around 0.5% or so for Im{uj). The second overtone showed much 
higher discrepancies, especially for Im{uj), with differences around 15% is some cases, while 
for Re{uj) this difference usually hovers around 1 — 2% (with two exceptions for £ = 8, when 
it was much bigger - almost 8% for the axial field). 
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n 


^NUM 


'^WKBl 












0.1105 -0.1008i 


0.1046 - 0.1152? 


1 





0.292939 - 0.0976663i 


0.2929 - 0.0977i 


0.2911 - 0.0980? 


1 


1 




0.264 - 0.307i 


0.2622 - 0.3074? 


2 





0.483644 - 0.0967590Z 


0.4836 - 0.0968Z 


0.4832 - 0.0968? 


2 


1 


0.470 - 0.300i 


0.4638 - 0.2956i 


0.4632 - 0.2958i 


2 


2 




0.4317 - 0.5034i 


0.4317 - 0.5034? 


3 





0.675367 - 0.0964997i 


0.675366 - 0.0965006Z 




3 


1 


0.667 - 0.288i 


0.660671 - 0.292288? 




4 





0.867417 - 0.0963923i 


0.867416 - 0.0963919i 




4 


1 


0.859 - 0.287i 


0.855808 - 0.290877i 




5 





1.059614 - 0.0963337Z 


1.05961 - 0.0963368? 




5 


1 


1.050 - 0.283i 


1.05004 - 0.290154? 




6 





1.251891 - 0.0963060i 


1.25189 - 0.0963051? 




6 


1 


1.224 - 0.283i 


1.24375 - 0.289736i 




7 





1.444214 - 0.0962866Z 


1.44421 - 0.0962852? 




7 


1 


1.433 - 0.282? 


1.43714 - 0.289473? 




8 





1.636565 - 0.0962724Z 


1.63656 - 0.0962719? 




8 


1 


1.629 - 0.288i 


1.63031 - 0.289297i 




8 


2 


1.584 - 0.422i 


1.61797 - 0.483757? 




9 





1.828939 - 0.0962639Z 


1.82893 - 0.0962626? 




9 


1 


1.822 - 0.288i 


1.82333 - 0.289173?? 




9 


2 


1.791 - 0.404i 


1.81225 - 0.483235? 




10 





2.021329 - 0.0962568i 


2.02132 - 0.0962558? 




10 


1 


2.014 - 0.288i 


2.01625 - 0.289083? 




10 


2 


1.990 - 0.-105/ 


2.0062 - 0. 182851/ 




11 





2.213730 - 0.0962522Z 


2.21372 - 0.0962507? 




11 


1 


2.207 - 0.288i 


2.20909 - 0.289015i 




11 


2 


2.174 - 0.417i 


2.19989 - 0.482567? 




12 





2.406139 - 0.0962487i 


2.40613 - 0.0962467? 




12 


1 


2.400 - 0.288i 


2.40186 - 0.288963? 




12 


2 


2.371 - 0.440i 


5 

2.39338 - 0.482347? 







n 


I^NUM 


I^WKBl 




1 





0.248229 - 0.0924905i 


0.2482 - 0.0926? 


0.2459 - 0.0931? 


1 


1 




0.2143 - 0.2941? 


0.2113 - 0.2958? 


2 





0.457595 - 0.0950044Z 


0.4576 - 0.0950? 


0.4571 - 0.0951? 


2 


1 


0.440 - 0.290? 


0.4365 - 0.2907? 


0.4358 - 0.2910? 


2 


2 




0.4009 - 0.5017i 


0.4023 - 0.4959? 


3 





0.656899 - 0.0956165i 


0.6569 - 0.0956? 


0.6567 - 0.0951? 


3 


1 


0.648 - 0.286? 


0.6417 - 0.2897? 


0.6415 - 0.2898? 


3 


2 




0.6138 - 0.4921? 


0.6151 - 0.4901? 


3 


3 




0.5814 - 0.6955i 


0.5814 - 0.6955? 


4 





0.853096 - 0.0958605i 


0.853095 - 0.0958601? 




4 


1 


0.844 - 0.285? 


0.841267 - 0.289315? 




5 





1.047915 - 0.095982H 


1.04791 - 0.095981? 




5 


1 


1.039 - 0.282? 


1.03822 - 0.289104? 




6 





1.242000 - 0.0960523i 


1.242 - 0.0960512i 




6 


1 


1.232 - 0.282? 


1.23379 - 0.288982? 




7 





1.435647 - 0.0960959Z 


1.43564 - 0.0960947? 




7 


1 


1.424 - 0.281? 


1.42852 - 0.288906? 




8 





1.629012 - 0.0961250i 


1.629 - 0.0961236i 




8 


1 


1.621 - 0.288? 


1.62272 - 0.288855? 




8 


2 


1.577 - 0.425? 


1.61032 - 0.483028? 




9 





1.822180 - 0.0961452Z 


1.82217 - 0.0961439? 




9 


1 


1.815 - 0.288? 


1.81655 - 0.288819? 




9 


2 


1.788 - 0.407? 


1.80542 - 0.48265? 




10 





2.015214 - 0.0961596Z 


2.01521 - 0.0961587? 




10 


1 


2. GOG - 0.288/ 


2.01012 - 0.288793/ 




10 


2 


1.984 - 0.410? 


2.00003 - 0.482374? 




11 





2.208148 - 0.0961714i 


2.20814 - 0.0961697i 




11 


1 


2.199 - 0.285? 


2.20349 - 0.288773? 




11 


2 


2.161 - 0.442? 


2.19426 - 0.482167? 




12 





2.401004 - 0.096180? 


2.40099 - 0.0961782? 




12 


1 


2.395 - 0.288? 


6 

2.39672 - 0.288758? 




1 o 


o 











n 


I^NUM 






2 





0.37367 - 0.08896i 


0.373691 - 0.088891i 


0.3732 - 0.0892? 


2 


1 


0.352 - 0.272Z 


0.346297 - 0.27348Z 


0.3460 - 0.2749? 


2 


2 




0.2985 - 0.4776i 


0.29852 - 0.47756? 


3 





0.599444 - 0.092703H 


0.599443 - 0.0927025Z 


0.5993 - 0.0927? 


3 


1 


0.587 - 0.278i 


0.582642 - 0.28129i 


0.5824 - 0.2814? 


3 


2 




0.551594 - 0.479047Z 


0.5532 - 0.4767? 


4 





0.809180 - 0.0941643i 


0.809178 - 0.094164H 


0.8091 - 0.0942? 


4 


1 


0.797 - 0.279Z 


0.796631 - 0.284334Z 


0.7965 - 0.2844? 


4 


2 




0.7727 - 0.4799i 


0.772695 - 0.4799? 


5 





1.012297 - 0.0948713i 


1.0123 - 0.0948706i 




5 


1 


1.002 - 0.279Z 


1.00222 - 0.285817Z 




6 





1.212013 - 0.0952667i 


1.21201 - 0.0952659Z 




6 


1 


1.203 - 0.286Z 


1.20357 - 0.28665Z 




7 





1.409741 - 0.0955106i 


1.40974 - 0.0955096ii 




7 


1 


1.401 - 0.286Z 


1.40247 - 0.287164i 




7 


2 


1.361 - 0.432Z 


1.38818 - 0.480709Z 




8 





1.606202 - 0.0956724i 


1.60619 - 0.0956707Z 




8 


1 


1.594 - 0.280i 


1.59981 - 0.287504i 




8 


2 


1.474 - 0.515i 


1.58721 - 0.480804i 




9 





1.801801 - 0.095784i 


1.80179 - 0.0957828Z 




9 


1 


1.795 - 0.288Z 


1.7961 - 0.287741? 




9 


2 


1.770 - 0.440Z 


1.78483 - 0.48087i 




10 





1.996796 - 0.0958649i 


1.99679 - 0.0958639i 




10 


1 


1.990 - 0.286Z 


1.99165 - 0.287912Z 




10 


2 


1.9G7 - 0.408/ 


1.98115 - 0.-18087/ 




11 





2.191345 -0.0959263i 


2.19133 - 0.0959245? 




11 


1 


2.185 - 0.287i 


2.18665 - 0.28804i 




11 


2 


2.154 - 0.415i 


2.17734 - 0.480952i 




12 





2.385555 - 0.0959724i 


2.38554 - 0.095971? 




12 


1 


2.379 - 0.287Z 


2.38124 - 0.288138? 




12 


2 


2.353 - 0.446i 


7 

2.37268 - 0.480979? 





Sch BH, EM L=8 Field 



M= 1.0ni,h = 0.01m 




t(m) 

FIG. 1: Fundamental mode and first overtone, for the electromagnetic £ = 8 field. The second 
overtone is much clearer than in the previous figures. 

APPENDIX A: FINDING THE OVERTONES NUMERICALLY 

A comment is due on the method of extracting the secondary (and, whenever apphcable, 
the tertiary) mode: we took the first (or dominant) mode {n = 0) and apphed an oscillatory 
fitting to it, in order to extract the data on its oscillatory frequency lor, amplitude A and 
decay rate uj. These data were taken with the maximum number of significant digits possible 
(usually 8 or 9), and we subtracted the fitted function, which had the form 

^'o = Aexp{iuj%t) exp(-u;°t), (Al) 

so that we could see the remainder of this operation. In all cases (except for the i = 0,1 
scalar and i = 1 EM field), we could find a remainder of the form 

^R = B exp(icj]ji(;) exp(-cj]t) + A, (A2) 

in which uji ^ Wq, characterizing a secondary mode (the first overtone), and with similar 
amplitudes {A^ B). The term A is just the \l/o term, but with a much reduced amplitude, 
indicating that the fitting operation has its precision limitations. This can be seen in the 
form of parallel curve envelopes in the figures from above. 
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Sch BH, Scalar L=7 Field 

M= 1.0m, h = 0.01m 




t(m) 

FIG. 2: Fundamental mode and first overtone, for the scalar £ = 7 field. The second overtone also 
appears, though not in a very evident manner. 

A similar procedure was adopted to find the second overtone, this time applying the 
oscillatory fitting to the first overtone, in which a new remainder similar to that seen in 
(IA2p was seen. In this latter case, however, only the very high ^-values could yield any 
valuable result when an oscillatory fitting was applied to it. 

Now it is the time to estimate the precision of the first fitting. As we can see again in 
the aforementioned figures, there is a given time tc for which the secondary mode and the 
dominant mode are roughly equal in magnitude, that is 

B exp{iLj]^t) exp{—uj]t) ^ S{Aexp{iuj^t) exp{—ujjt)). (A3) 

The variation of the dominant mode can be expressed as 

S{Aexp(iuj^t) exp(—uJit)) = exp(iu;^t) exp(— cj/t) + 

+ Ai exp{iuRt) exp{-ujjt) [5uj% + i5(J^'^) 

In what follows, we assume 6A to be much smaller than A. Such approximation implies 

B eyip(iiu]^t) exp(— tujt) ~ Ait5to exp^iuj^) exp{—uj]t) , (A5) 

in which we have used 6uj = 5a;^ + i(5a;j. Working with the moduli of the quantities above, 

9 



we arrive at 



M^q r^' (A6) 

Upon substituting t = tc and noting that 1^41 « \B\, we arrive at 

Such an estimate was done for each field and for each £- value, and \6lj\ was compared to |a;|, 
so as to determine the approximate degree of precision in the computation of the dominant 
mode frequency. A similar estimate can also be done, in principle, for the first overtone 
fitting precision. We have determined the ratio m — \^\ - for the n — mode - to be 
smaller than 10~^, in all cases we have dealt with, except for the axial £ — 2 case, for which 
m w 10~^ (hence its smaller number of significant figures). 
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